Interpreting the Behavior of Time Resolved Gamma-Ray Burst 

Spectra. 
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ABSTRACT 

In this paper, we explore time resolved Gamma-Ray Burst (GRB) spectra in the context 
of the synchrotron emission model presented in Lloyd and Petrosian (2000; LP00). First, we 
show that our model - which involves three distinct emission regimes - can provide excellent fits 
to the time resolved spectra of GRBs, and we present these results for a few bursts. We then 
describe how the phenomenological Band spectrum (Band et al., 1993) can be interpreted in 
the context of our models based on the value of the low energy photon index a. We discuss the 
types of correlations one would expect to observe among the Band parameters if these models 
are correct. We then compare these predictions to the existing data, combining a sample of 
2,026 time resolved spectra (from approximately 80 bursts). We show that the correlations 
found in the data are consistent with the models, and discuss the constraints they place on 
the emission physics. In particular, we find a (~ ia) negative correlation between the peak of 
the vF v spectrum, E p , and the low energy photon index a for bursts with —2/3 < a < 0, in 
contrast to what is predicted by the instrumental effect discussed in LP00. We suggest that 
this correlation is simply due to the mechanism responsible for producing en's above the value 
of —2/3 - namely, a decreasing mean pitch angle of the electrons. We also show that E p is 
correlated with the photon flux, and interpret this as a result of changing magnetic field or 
characteristic electron energy between emission episodes. Finally, we discuss the implications 
our results have on particle acceleration in GRBs, and prospects for further testing these models 
with the anticipated data from HETE-2, Swift and GLAST. 



1. Introduction 

Except for a few isolated bursts (see, for 
example, Tavani, 1996, Brainerd, 1996), most 
of the analyses of the prompt spectral data of 
Gamma-Ray Bursts have employed phenomeno- 
logical models - in most cases, using the so-called 
Band spectrum (Band et al., 1993). This model 
is essentially a smoothly broken power with a 
low energy photon spectral index a, a high en- 
ergy photon index /?, and a break energy E p . 
It has been shown that such a model describes 
most GRB spectra very well (e.g. Mallozzi ct 
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al., 1996, Lloyd & Petrosian, 1999, Preece et al., 
1999). There have been some attempts to ex- 
plain or interpret some of the global properties 
of these spectral parameters in terms of a physi- 
cal model (e.g., Totani, 1999, Preece et al., 1998, 
Ghissilini et al., 2000), with inconclusive results. 
In Lloyd and Petrosian (2000; herafter LP00), 
we showed that realistic synchrotron models can 
qualitatively explain the global distributions of 
the time averaged Band spectral parameters. Our 
model modified the usual simple picture of opti- 
cally thin synchrotron emission from a power law 
distribution of electrons with a sharp low energy 
cutoff, by accounting for: 1) the possibility of a 
smooth cutoff to the low energy electron distri- 
bution, 2) radiation from an anisotropic electron 
distribution with a small mean pitch angle, 3) 
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synchrotron self-absorption, and 4) the impor- 
tant instrumental effect in which the value of the 
fitted parameter a decreases as E p approaches 
the lower edge of the BATSE window. We have 
envisioned a realistic scenario in which particle 
acceleration and synchrotron losses occur continu- 
ally and simultaneously behind an internal shock 
(which produces a single pulse in the GRB time 
profile) with the characteristic acceleration time 
shorter than the loss time, so that synchrotron loss 
effects are only evident in the particle distribution 
spectrum at energies much larger than what is 
relevant for our discussion here (at energies where 
the inequality is reversed and the loss time be- 
comes shorter than the acceleration time). This 
model stands up well to the global distributions of 
GRB parameters. In particular, it can accommo- 
date the bursts with a above the so-called "line 
of death" value a = —2/3 (Preece et al., 1998). 
However, the tests conducted so far have involved 
time averaged properties of the bursts. Since each 
pulse is regarded as a separate emission episode 
(or internal shock) in our model, the averaging 
over many pulses may obscure the real physics 
of each episode. It would, therefore, be useful to 
examine the time resolved spectral properties of 
GRBs, so that one may compare the values and 
the correlations between the spectral parameters 
from pulse to pulse and perhaps within a pulse, 
in hopes of gaining insight on the evolution of 
physical parameters throughout a GRB. 

The temporal evolution of GRB spectral pa- 
rameters has been studied by several authors (e.g. 
Norris ct al., 1986, Kargatis et al., 1994, Ford et 
al., 1995, Crider et al., 1997, and Preece et al., 
1998, Ryde & Svensson, 2000). These studies have 
attempted to look for global trends in the data - 
in particular, Norris et al. and Ford et al. re- 
port a hard-to-soft evolution trend, while Crider 
et al. report both a hard-to-soft and "tracking" 
trend (in which one or more spectral parameters 
track the time profile of the burst). The tempo- 
ral behavior has not been rigorously interpreted 
in terms of an emission scenario (although Crider 
et al. attempt to explain at least the hard-to-soft 
evolution as evidence of a Comptonized spectrum 
in which the medium is expanding, causing E p to 
decrease with time). Furthermore, with the excep- 
tion of Kargatis et al., these studies have perhaps 
over-simplified the extremely varied and complex 



evolutionary behavior seen in most GRBs. 

The purpose of this paper is to use our syn- 
chrotron emission models as diagnostics for inter- 
preting the time resolved GRB spectral parame- 
ters. Using a physical model - rather than a phe- 
nomenological model such as the Band spectrum 
- to characterize the GRB spectra, allows us to 
gain insight into the evolution of actual physical 
parameters in the GRB plasmas. For example, in 
the case of optically thin emission by an isotropic 
distribution of electrons, E p cx ^B±, where j m is 
the characteristic cutoff or turnover energy of the 
electron distribution, and B± is the perpendicu- 
lar component of the magnetic field. The evolu- 
tion of E p reflects the evolution of 7™ and/or B±. 
Clearly, just examining the evolution of this one 
spectral parameter will not lead to insights on the 
evolution of a single physical quantity, since the 
former (E p ) depends on more than one physical 
variable (j m , B±). However, examining the evo- 
lution of this parameter and its correlation with 
other spectral parameters can help break some of 
the degeneracies and the temporal behavior of the 
plasma parameters can begin to be elucidated. 

Furthermore, once we get a handle on these pa- 
rameters and their evolution throughout the GRB 
event, we can begin to gain some insight into the 
fundamental problem of particle acceleration in 
GRBs. The particle acceleration and loss (e.g. to 
to radiation and otherwise) mechanisms determine 
the emitting electron distribution. The physics be- 
hind the particle acceleration determines the value 
of the minimum and maximum electron Lorentz 
factors, the smoothness of the low energy cutoff to 
the electron distribution, the high energy power 
law index p of the electron distribution, and is 
of course intimately related to the magnetic field. 
The presence of a significant number of bursts with 
steep high energy spectra - that is, with a high en- 
ergy photon index (3 = — (p+l)/2 less than about 
—3 - suggests that p can be significantly larger 
than the so-called "universal index" p w 2.2 due 
to first order Fermi acceleration in a relativistic 
shock (see, e.g., Kirk et al., 2000 and references 
therein); thus, a different type of particle accelera- 
tion mechanism may be at work. Furthermore, the 
evidence for the presence of synchrotron radiation 
from electrons with small pitch angles (see LP00 
and below) also suggests that the usual assump- 
tion of isotropicization of the particle spectrum on 
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short timescales may not be warranted. It is clear, 
then, that a burst's spectrum can help guide us in 
a detailed study of particle acceleration in internal 
shocks. 

In this paper we will further develop and ex- 
plain our synchrotron emission models, and then 
use them to examine the behavior of the time re- 
solved GRB spectra in the context of an in internal 
shock scenario. In §2, we review the characteristics 
of our synchrotron models, and show that these 
models provide good fits to the existing data. In 
§3, we discuss the types of correlations expected 
among the spectral parameters, in the context of 
our models. In §4, we examine correlations present 
between spectral parameters from Band fits car- 
ried out by Preece et al. (1999) in a sample of 
2,026 time resolved spectra. We find a strong pos- 
itive correlation between the total photon flux,/ 7 , 
and the peak of the vF v spectrum, E p , and cor- 
relations between E p and the low energy photon 
index, a, which differ depending on whether a is 
above or below the value of ~ —2/3. We inter- 
pret these results, in the context of our models, 
as reflecting changes in the physical parameters 
from one emission episode to the next (one inter- 
nal shock to the next). In §5, we present a few 
cases of spectral evolution for individual bursts, 
simply to illustrate how the models can be used 
to infer something about the physics in a partic- 
ular burst (from shock episode to shock episode). 
Finally, in §6, we present conclusions and discuss 
the implications of our results on particle acceler- 
ation in GRBs. 

2. Some Distinct Synchrotron Emission 
Scenarios 

The details of our synchrotron emission models 
are described in LP00, where it is shown that the 
low energy spectral index a plays a key role in un- 
derstanding the emission mcchanism(s) at hand. 
In this section, we describe three possible emis- 
sion scenarios with distinct asymptotic low energy 
spectral behavior. 

2.1. Isotropic Pitch Angle Distribution of 
Electrons (IPD) 

This is the familiar optically thin synchrotron 
emission from a power law electron energy spec- 
trum, with an isotropic pitch angle distribution; 



but, in contrast to most analyses, here we con- 
sider an electron distribution with a smooth low 
energy cutoff: AT (7) oc pp^^ra • Note that 
for high energies (7 > -f m ), the spectrum goes as 
7 _p , while for low energies (7 < 7 m ), the spec- 
trum goes as 7 9 . Hence, q denotes the steepness 
of the electron low energy cutoff (note that an 
actual "cutoff", in the sense that AT (7) — ► as 
7^0, requires q > 0). The asymptotic behavior 
of the synchrotron (photon number) spectrum for 
q > -1/3 is: 

p - J "~ 2/3 V « v m = |^ssin*7^ , , 

where F 1 is the photon flux, ^ is the electron pitch 
angle, and vb = where B is the magnetic field. 
Note that the peak of the vF v spectrum will occur 
at E p oc v. m oc _Bsin\D'7^ l , and that the aymptotic 
low energy index below this break is a — —2/3. 
We point out that if q < —1/3, the low energy 
asymptotic index is a = (q — l)/2; since q < —1/3 
does not constitute a "cutoff", we do not discuss 
this further and limit our discussion to cases with 
< q < 00. However, even though the low energy 
asymptotic index of the photon spectrum is always 
—2/3 for these (latter) cases, this does not mean 
that the value of q does not play an important role 
in the observed spectrum. As shown in LP00 and 
discussed in §3, the smaller the value of q is, the 
lower the frequency at which the asymptotic value 
of —2/3 is reached; because of the finite width of 
the detector spectral window, this can cause the 
fitted value of a to be significantly less than —2/3. 

2.2. Small Pitch Angle Distribution of 
Electrons (SPD) 

This spectrum results from optically thin syn- 
chrotron emission by electrons with a mean pitch 
angle >t < 1; the analysis of synchrotron radiation 
in this regime was first done by Epstein (1973). 
For high density, low magnetic field plasmas, the 
Alfven phase velocity is less than the speed of light 
and (therefore) the speed of the relativistic parti- 
cles under consideration here. In this case, the 
pitch angle diffusion rate of the electrons interact- 
ing with plasma turbulence is much larger than 
the acceleration rate; consequently, the acceler- 
ated electrons will have an isotropic pitch angle 
distribution. However, for the low density, high 
magnetic field conditions expected for the sources 
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of GRBs, the opposite is true. In this case the 
amplitude of the electric field fluctuations exceeds 
that of the magnetic field so that the above sit- 
uation is reversed (see e.g. Dung and Petrosian, 
1994 and Pyradko and Petrosian, 1998). Then 
the pitch angle distribution of the accelerated elec- 
trons could become highly anisotropic as required 
in the small pitch angle model. The shape of this 
spectrum depends on just how small the pitch an- 
gle is. For ^ <C 1, but ^7 TO ~ 1, we have: 

( v>° v«.v s = |^ B /( 7m * 2 ) 

F-r = { v s « v « v m (2) 

There are two breaks in this spectrum - one at 
v, m and one at v s . Because the Band spectrum 
can only accommodate one break, spectral fits to 
this model will put the parameter E p at one or 
the other of these two breaks, but most likely at 
v m because for p > 5/3 (or for high energy pho- 
ton index (3 < —4/3 which is the case for most 
bursts), the break across v m is more pronounced 
than across v s . In this case, the low energy photon 
index a will fall somewhere between —2/3 and 0. 

However, as the pitch angle \& decreases such 
that * <C l/7mi then the v~ 2 l 3 portion of the 
spectrum disappears, and only the t>° portion is 
left. In this case we have: 

Tp _ J f° K<f s = §f£!7m ,o\ 

where E p oc Bj m (see Epstein, 1973 for a more de- 
tailed description of the behavior of the spectrum 
in this regime). [We note that Medvedev (2000) 
has developed a model in which the transverse de- 
flections of electrons in highly non-uniform, small 
scale magnetic fields are smaller than the elec- 
trons' relativistic beaming angles (<~ l/7 e ), so that 
the entire trajectory of the electron is observable. 
In this case, a so-called "jitter" spectrum is ob- 
tained (Medvedev, 2000) which has some of the 
same low energy characteristics as the SPD spec- 
trum - in particular, the low energy photon index 
in this model also has a value of 0.] 

2.3. Self-Absorbed Spectrum (SAS) 

If the magnetic field and density are such that 
the medium becomes optically thick to the syn- 
chrotron photons with frequency v < u a , then, for 



v a <v m , we have the following spectrum: 

{v 1 v < v a , 

v- 2 ^ v a <^v<^v m ,. (4) 

z,-(p+D/ 2 v m ^>v 

In that case, E p oc v a ~ 10(nZ) 3/5 5 2/5 7m 8/5 r 9/5 
Hz, where I and n are the path length and particle 
density in the co-moving frame, and we have as- 
sumed an electron energy distribution index p = 2. 
[We have also assumed an isotropic distribution 
of electron pitch angles; for a small pitch an- 
gle distribution, the i/~ 2 / 3 portion in equation 
(4) would be replaced by v° .] For v a > v m we 
just have one break at v a with a low energy pho- 
ton index of a = 3/2 (in both the isotropic and 
small pitch angle cases). The possibility of self- 
absorption in GRBs is a controversial issue. We 
have shown (LP00) that there are bursts for which 
a self-absorbed spectrum is a better fit than an 
optically thin one. We also found that in these 
cases, the absorption frequency tends to be near 
the lower edge of the BATSE window. In addi- 
tion to this, Strohmeyer et al. (1998) found that 
a number of bursts observed by GINGA with E p 's 
in the range 2 to 100 keV have steep (a <~ 1) 
low energy spectral indices consistent with a self- 
absorbed spectrum. This raises interesting ques- 
tions about the physics of the ambient plasma, 
because self-absorption in a GRB requires fairly 
large magnetic fields and particle densities. For 
example, if the absorption frequency is less than 
the minimum electron frequency, the optical depth 
to synchrotron self-absorption is 

t - (Z/10 13 cTO)(n/10 8 cm- 3 )(B/10 8 G) 2/3 

(7m/50)- 8/3 (r/10 3 ) 3 (^ ofcs /40fcey)- 5 / 3 (l+z)- 5 / 3 , 

(5) 

where u (, s is the absorption frequency in the ob- 
server's frame, and z is the redshift of the GRB. 
Note that this frequency falls within BATSE's 
spectral window under certain, perhaps somewhat 
extreme, conditions. The physical processes re- 
quired to achieve these conditions will need to 
be theoretically established if the data prove self- 
absorption to be a viable model. We point out 
that the next generation of GRB dedicated tele- 
scopes - namely SWIFT and HETE-2 - will obtain 
more spectral data in energy ranges lower than the 
BATSE threshold (of about 25 keV) and can firmly 
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establish the presence or absence of a self-absorbed 
portion of the low energy spectrum. 

Figure 1 shows the various spectra for the dif- 
ferent emission regimes. Note that we have plot- 
ted the F v spectrum rather than F 1 — F v /u to 
emphasize the differences in the various low en- 
ergy slopes of the spectra between the different 
emission regimes. Throughout the rest of the pa- 
per, we use the value of the low energy photon 
index a to distinguish between the different emis- 
sion scenarios, where the IPD case is defined 
by a < -2/3, SPD by -2/3 < a < 0, and SAS 
by a > 0. 




hy(keV) 

Fig. 1. — Various synchrotron energy spectra, 
F v (in arbitrary units), as a function of energy 
hv in kcV. The dot-dashed line is optically thin 
radiation from an isotropic distribution of elec- 
trons (IPD) with a sharp minimum energy cut- 
off, while the dotted and short-dashed lines show 
the IPD synchrotron spectra for smooth cutoffs 
to the electron distribution (q = and q = 2 re- 
spectively). Note that for finite values of q, the 
asymptotic v 1 ^ spectrum is achieved at photon 
energies much lower than E p . The solid line shows 
a self-absorbed spectrum (SAS) for v a > v m . 
The long-dashed and dot-dashed lines indicate the 
small pitch angle distribution (SPD) case for small 
($ C 1) and very small (7* <C 1) pitch angles, 
respectively. The vertical lines mark the approxi- 
mate width of the BATSE spectral window. 



2.4. Spectral Fits - Directly Testing the 
Model 

To test the how well these synchrotron mod- 
els actually describe the existing data, we have 
performed time resolved spectral fits to a sample 
of data from the BATSE archive. We use 128 
channel, 128ms time resolved HER data, which 
is obtained for the most brightly illuminated of 
the eight detectors from the on-line archive at: 
cossc.gsfc.nasa.gov/compton /data/batse/trigger. 
We plot up the total counts as a function of time 
(the burst time profile) summing over all energy 
bins; from this, we pick out time intervals for the 
background and over which to do our spectral fits. 
We then subtract off the background counts aver- 
aged over our specified "background" time inter- 
vals from our raw spectral counts data for each 
energy bin (the alternative method of subtract- 
ing a fitted background photon model off of the 
spectral photon model is mathematically equiv- 
alent - that is, it is equivalent to subtract off 
the background before convolving with the detec- 
tor response matrix (DRM) or after convolving 
with the DRM ). We then fit our photon models 
to the data by convolving them with the DRM 
to get model counts and then minimized x 2 = 
(data-model counts) 2 /er 2 , via a downhill simplex 
method. 

Our models are described in detail in LPOO and 
briefly above. We fit each spectrum to all 3 emis- 
sion scenarios and then evaluate the fits based on 
their values of a reduced x 2 - In Figure 2, we show 
examples of spectral fits in each emission regime. 
Each fit is taken at a time during the burst spec- 
tral evolution when the a parameter corresponded 
to the respective model. For example, in the top 
panel - burst 1663 - the spectrum is from a time 
when a ~ —2/3, while in the middle panel - burst 
105 - the spectrum is from a time in the profile 
when a = 0. Similarly, for the bottom panel, this 
spectrum corresponds to a time when a = 1. The 
reduced x 2 are 0.34, 0.33, and 0.50 for the top, 
middle and bottom panels respectively. In gen- 
eral the best model turns out to correspond to 
the emission regime suggested by Band's a values, 
which confirms our proposed method of physically 
interprctating Band fits based on the bursts' low 
energy photon index (for example, an IPD fit to 
the spectrum of burst 105 gave a x 2 > 1 compared 
to the x 2 = 0.33 for an SPD fit). However, we 
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Fig. 2. — Fits to the three different synchrotron 
models described in §2. The emission regime that 
best fits each burst (indicated in each figure) is 
the regime inferred from the Band fit value of a 
at the time of each spectrum. [Note that these are 
the actual counts spectra of the burst, which is 
the photon spectrum convolved with the detector 
response matrix.] 

have not yet been able to confirm this for the self- 
absorbed case; although in some bursts the SAS 
was clearly better than the IPD case, SPD seemed 
to do statistically as well for the small sample of 
bursts we have tried. As mentioned above, better 
low energy data is needed to definitively establish 
the existence of a self-absorbed component in some 
GRB spectra. Nonetheless, in this paper, we inter- 
pret those bursts with a low energy photon index 
larger than the SPD limit of a = as SAS cases. 

3. Expected Correlations Among Spectral 
Parameters 

Ideally we would like to carry out such fits 
to a large sample of bursts, and characterize the 
physics of each burst as well as trends among large 
samples of bursts directly through the physical pa- 
rameters yielded by the fits. However, this is a 
large and time-consuming task, and in fact we can 
learn a great deal about the underlying physics 
of GRBs through an analysis of the phenomeno- 



logical spectral parameters in terms of a physi- 
cal model. We have already shown that the Band 
parameter a is a good indication of the relevant 
synchrotron emission regime (IPD, SPD, or SAS). 
Examining correlations among the various Band 
parameters can lead to additional insights on the 
physics governing the emission, when interpreted 
in the context of a physical model. Below, we 
discuss the types of correlations one might expect 
among the Band spectral parameters for the dif- 
ferent emission scenarios described above. In §4 
and 5, we compare these correlations with what 
we find in the GRB data. 

3.1. a — E p Correlation 

Two different correlations are expected for 
these parameters. 

1. We expect a positive correlation between E p 
and a due to the instrumental effect described in 
LP00. If E p is close to the edge of the BATSE 
window, the low energy photon index may not yet 
have reached its asymptotic value and a smaller 
(or softer) value of a (relative to the asymptotic 
value) will be determined. A smooth cutoff to the 
electron energy distribution will exacerbate this 
effect because for a smoother cutoff (or a lower q), 
the low energy asymptote is reached farther away 
from E p , nearer to (or even below) the low en- 
ergy edge of the detector spectral window. Note 
that a dispersion in the smoothness of the low en- 
ergy cutoff will tend to wash this correlation out to 
some degree, as seen in Figure 4 of LP00. For the 
cases of small pitch angle radiation and the self- 
absorbed spectrum, this effect will be weaker be- 
cause the low energy asymptotes are reached more 
quickly (i.e. at energies closer to E p ) than for the 
isotropic optically thin case (see Figure 1). 

2. We also expect evidence of a negative corre- 
lation between E p and a as we transition from the 
IPD to the SPD regime, i.e. for -2/3 < a < 0. 
In this case, the pitch angle decreases so that 
E p oc sin^P decreases, if all other physical parame- 
ters (B and j m ) remain constant. In addition, as 
we go from the small pitch angle regime, \D r 7 m ~ 1 
("J <C 1), to the very small pitch angle regime, 
to ^7 TO <C 1, the z^~ 2 / 3 portion of the spectrum 
disappears, and we are left with only the v° por- 
tion. In other words, as the mean of the pitch 
angle distribution decreases to very small values, 
E v decreases and the value of a decreases from 
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—2/3 to 0. This negative correlation will compete 
with the positive instrumental correlation men- 
tioned above. 

3.2. 0-Ep Correlation 

We expect a similar correlation between E p and 
(3 due to instrumental effects. A dispersion in the 
high energy electron index p will tend to reduce 
this correlation. However, in practice - partly 
because the high energy spectral data are not 
very constraining and partly because E p is usually 
well below the upper edge of the BATSE window 
(about 1.5 MeV) for those bursts with spectral fits 
- this correlation between E p and (3 is not evident 
in the data. For a few sample spectra, we find that 
E p has to be greater than around 1100 keV before 
[3 is affected by this instrumental effect. 

3.3. Total flux - E p Correlation 

We might expect a positive correlation between 
E p and the flux of the burst. If E p changes ei- 
ther due to a change in the magnetic field or 7 TO , 
then the flux, which also depends on positive pow- 
ers of both of these parameters, will also increase 
(Pacholczyk, 1970). This effect of course will be 
weakened to some degree by the distribution of 
redshifts of GRBs (if we examine the whole spec- 
tral sample instead of one burst). However, as 
shown in Lloyd et al., 2000, the cosmological con- 
tribution to such a correlation (higher redshift re- 
duces the observed value of E p and flux for a given 
burst) is negligible due to the large intrinsic dis- 
persion in the luminosity function and intrinsic E p 
distribution, so that any correlation we do observe 
can be attributed to an intrinsic effect. 

4. The Time Resolved Spectral Data - 
Global Behavior 

We showed in LP00 that the time averaged 
spectral parameter distributions arc consistent 
with the models of synchrotron emission described 
above. We also showed for several GRBs, that 
the evolution of different spectral parameters of a 
particular burst track eachother throughout their 
time evolution in a way that is easily interpreted 
in our models. We want to test our models and 
in particular investigate the latter point in more 
detail by examining the behavior of a large sample 
of time resolved spectral parameters. The ideal is 



to learn something about how the plasma param- 
eters are changing between emission episodes in 
a burst, by examining the time resolved of the 
Band spectral parameters in the context of our 
synchrotron models. As mentioned in §3, because 
the spectral parameters can depend on more than 
one physical quantity, we look for correlations 
between the spectral parameters (rather than ex- 
amining the evolution of one spectral parameter 
in time), in order to break some of the degenera- 
cies in interpreting the evolution of the physical 
parameters. 

One way to do this is to look for particular 
trends between pulses within individual bursts. 
For example, Crider et al. (1997) claim to see 
a positive correlation between E p and a in 47 
individual bursts 3 . The advantage of looking at 
spectral evolution within individual bursts is that 
correlations between the spectra do not contain 
any dispersions or contributions that might arise 
due to redshift effects. However, looking for global 
evolutionary trends in the data by examining in- 
dividual bursts is a difficult and challenging task, 
particularly when we differentiate between differ- 
ent emission regimes. This is not only because 
there are a small number (typically ~ 20) of time 
resolved spectra per burst, but also because there 
are a small number of points per emission regime, 
particularly for the SPD and SAS regimes. There- 
fore, attaching a significance to correlations be- 
tween spectral parameters in different emission 
regimes for a single burst is in general not sta- 
tistically robust (Efron, private communication) 
and there is therefore no reliable way to compare 
results with other bursts in order to establish gen- 
eral trends in the evolution. 

We would like to investigate if there are any 
average trends present among the spectral param- 
eters in each emission regime. Although we do 
examine the behavior of some individual bursts in 
the next section, this goal is best accomplished by 
combining all 2,026 time resolved spectra avail- 
able, and searching for any global trends in this 

3 We note that they do not distinguish between different 
emission regimes in their analysis. They have also found ev- 
idence for a negative correlation in some individual bursts, 
but with admittedly low statistical significance. They do 
not report the a values for these bursts (so that they might 
be interpreted in terms of a particular "emission regime" 
of our models). 
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sample. Trends present in this entire sample will 
reflect the trends of temporal behavior in indi- 
vidual bursts on average, and can therefore help 
us gain insight on the evolutionary trends present 
in individual bursts. We can then test whether 
these trends are consistent with our models of syn- 
chrotron emission, and - if so - ideally learn some- 
thing about the evolution of the physical condi- 
tions (such as the magnetic field and Lorentz fac- 
tors) throughout a GRB. 

4.1. Data and The a Distribution 

Our data is taken from the catalog of Preece, et 
al. (1999), which contains high energy resolution, 
time resolved spectral fits to a large number of 
BATSE bursts (see their paper for discussion of 
data type, time and energy resolution, etc.). Our 
sample consists of individual spectra for which the 
HER data type is used in the spectral fit (because 
of its superior energy resolution; see Preece et al., 
2000 for description of data types), and for which 
the Band function provided a reliable fit in the ~ 
20 keV to ~ 1.5 MeV range. These criteria leave us 
with 2, 051 spectra. We then elimate 25 additional 
individual spectra, because the error bars on the 
spectral parameters are 0, indicating an error in 
the fitting procedure. This leaves us with 2,026 
spectra; although this is only a fraction of the ~ 
5, 000 spectra in the Preece catalog, we believe it is 
an accurate representation of at least those bursts 
which are described by the Band spectrum. And 
as discussed in Preece et al. (1999), this spectral 
form describes the large majority of bursts very 
well (see also Band et al., 1993, Mallozzi et al., 
1996, Lloyd and Petrosian, 1999) 4 . 



4 For some bursts, however, this is not the case. For exam- 
ple, Preece et al. (1996) showed that in a sample of about 
90 time averaged spectra, about 14% showed an X-ray ex- 
cess in the 7 — 20 keV range. Strohmeyer et al. (1998) also 
found an X-ray excess in at least one GINGA burst. Such 
an excess would tend to affect the low energy photon index 
- giving a softer (lower) value of a than if there were no ex- 
cess. This, in turn, would add some scatter to the expected 
correlations between a and E v discussed in §3.1. However, 
because we arc looking at the subset of data for which the 
Band spectrum provides a reliable fit, and we are examin- 
ing time resolved spectra between 20 — 1500 keV (above the 
range in which theses excesses have been found) , we do not 
expect this would significantly affect our results. Of course, 
one can always add additional emission components to any 
model; but because we have no evidence of it in our data, 
we take the simplest interpretation of a single component 
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Fig. 3. — Histogram of low energy photon index a, 
a compilation of 2,026 time resolved spectral fits 
from Preece et al., 1999. The different emission 
regimes and percent of spectra in each regime are 
marked on the figure. 

4-1.1. The a Distribution 

As discussed in §2, the parameter a is the 
best parameter for distinguishing between the var- 
ious synchrotron regimes. Figure 3 shows a his- 
togram of a (taken from the 2,026 time resolved 
fits of Preece et al., 1999) with each regime clearly 
marked 5 . For our sample, we find the percent of 
spectra in each regime is as shown in Table 1. Note 
that there are a significant number of spectra in 



emission model and see what we can learn from the data 
under this assumption. 
5 We point out that the peak of this distribution falls just 
slightly higher than the a distribution presented in Preece 
et al., 2000, which peaked at around — 1. This may be 
partly due to the differences in the way wc binned our data, 
but primarily reflects the fact that we have only included 
those bursts with Band fits in our sample; for example, 
bursts that are described well by a "Comptonized" spec- 
trum or simple broken power-law spectra tend to have a 
slightly lower value of a than those described by the Band 
spectrum, and may shift the peak of the a distribution to- 
ward lower values. In any case, we note the difference is 
within the average la error on a (< Aa >~ 0.28), and 
remind the reader we are only examining those bursts for 
which the Band function provided an acceptable spectral 
fit. 
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Regime 


a range 


# 


% 


IPD 


a < -2/3 


1122 


55 % 


SPD 


-2/3 < a < 


805 


40 % 


SAS 


< a < 3/2 


99 


5 % 



Table 1: Number spectra in each emission regime. 



the SPD regime and SAS regime. It is important, 
however, to briefly discuss how the error bars on 
the parameter a affect the interpretation of this 
distribution. First, we point out that Preece et al. 
(1999) showed that the error bar on a cannot alone 
account for the large dispersion in the distribution. 
In our sample, wc find that the average la error 
on a is < Aa >~ 0.28, suggesting (in agreement 
with Preece et al.'s more detailed analysis) that 
there is still a significant number of bursts above 
the line of death a = —2/3. To quantify this, we 
have computed the number of bursts in each emis- 
sion regime, using a's upper and lower la limits. 
If we take all values of a at their upper limits 
(a + Aa), wc find 42%, 46%, and 12% of bursts in 
the IPD, SPD, and SAS regimes respectively. Tak- 
ing all values of a at their la lower limits (a— Aa), 
we find 69%, 29%, and 2% of bursts in the IPD, 
SPD, and SAS regimes respectively. Although the 
error bars on a can make some difference as to 
the numbers of spectra in each regime, we will see 
that this does not affect the qualitative nature of 
our conclusions below. We now discuss the corre- 
lations present in the data and their consistency 
with what we expect in the context of the three 
synchrotron emission scenarios. 



data respectively (we do point out that the scat- 
ter of E p in each bin is not necessarily Gaussian; 
see, e.g., Preece, et al., 1996 and Preece, et al., 
1998b). The observed trends are consistent with 
what is expected from our model in each emission 
scenario, and tell us something important about 
the role various effects play in the correlations, as 
wc discuss below. 




-1.5 -1 -0.5 0.5 



-1 -0.5 0.5 



4.2. Observed a — E p Correlation 

Figures 4a and b show the binned average cor- 
relation between E p and a present in the time re- 
solved data. For each of these figures, we have 
sorted a in ascending order and binned the data 
every 100 points (the horizontal error bars indi- 
cate the size of the bins). We then computed the 
average and median E p for these 100 points. In 
Figure 4a, we show the average E p vs. a, where 
the vertical error bars arc simply the variance of 
the mean value of E p ; in Figure 4b, we have plot- 
ted the median E p , where the solid and dotted 
vertical error bars indicate the range of E p about 
the median value that includes 68% and 90% of the 



Fig. 4. — Peak of vF v spectrum E p vs. low energy 
index a for binned 2,026 time resolved spectral 
fits. The upper panel plots the variance of the av- 
erage E p as the vertical error bars, while the lower 
panel plots the median values of E p and the range 
that contains 68% (solid lines) and 90% (dotted 
lines) of the data around the median. Each regime 
of emission is marked on the plot. The dashed 
line in the SPD regime is an example of how E p 
changes as a function of a within the BATSE spec- 
tral window as the mean of the electron pitch angle 
decreases, but all other parameters (B, -f m , p, q) 
remain constant. 
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4.2.1. IPD Regime 

Performing a Kcndell's r test on all of the (un- 
binned) data, we find a 9a positive correlation be- 
tween a and E p in the IPD regime. To account 
for both the error in a and E p , we have performed 
this test on all permutations of correlations be- 
tween the lower and upper values of a (from the 
lcr error bars) with the lower and upper values of 
E p . In addition we have averaged the value of t 
from 100 sets of data, in which - for each data 
point - a and E p are drawn from Gaussian dis- 
tributions with means equal to the parameter val- 
ues given in the catalog and standard deviations 
corresponding to the error bars, for each of the 
data points, correlating these deviates, and tak- 
ing the average value for the correlation. In all 
cases, we find a highly significant (> 6a) correla- 
tion. The positive correlation between a and E p in 
the IPD regime can be simply understood by the 
instrumental effect discussed in §3 above. How- 
ever, for a given value of the low energy electron 
distribution cutoff parameter q, the correlation is 
expected to be much stronger than observed (see 
Figure 4 of LP00). Of course, any dispersion in 
q will tend to weaken and flatten this correlation. 
In fact, we can give a rough but quantitative es- 
timate of the q distribution required to produce 
the observed E p — a correlation seen in Figure 4. 
For each burst (under the assumption that opti- 
cally thin synchrotron emission from an isotropic 
distribution of electrons produces its spectrum), 
we can estimate the q value necessary to produce 
the burst's E p and a values, simply from the de- 
termined relationships between a and E p for dif- 
ferent values of q in Figure 4 of LP00. In Figure 
5, we present an estimate of the distribution of q 
required to reproduce the data and therefore the 
shallow correlation between a and E p observed in 
the IPD regime. Note that the dotted part of the 
histogram is not at all well constrained because 
the correlation dramatically weakens for very high 
values of q and the solutions become quite degen- 
erate (that is, for a given E p , a q of 5 or 10 may 
produce the same value of a). 

4.2.2. SPD Regime 

The sign of the correlation between a and E p re- 
verses for values of a > —0.7, which is suggestively 
very close to where we expect a transition from the 



Fig. 5. — Estimated distribution of q, given the 
observed correlation between a and E p in the IPD 
regime. There is a degeneracy among high values 
of q (that is, for different q 7 s > 5, we obtain the 
same observed a value for a given E p ); the dotted 
portion of the histogram attempts to account for 
this degeneracy and is simply a reasonable guess 
at how this portion of the histogram behaves. 



IPD to SPD regime. Performing a Kcndell's r test 
on all of the (unbinned) data in this regime, we 
find a Aa negative correlation between a and E p . 
Again, to account for the error in both a and E Pl 
we have performed this test on all permutations 
of correlations between the lower and upper val- 
ues of a (from the la error bars) with the lower 
and upper values of E p , as well as averaged the 
r value from 100 sets of data drawn from distri- 
butions based on the existing data, according to 
the prescription described in §4.2.1. In all of these 
cases, we find a significant (> 3a) negative corre- 
lation. As mentioned in §3, this type of correlation 
is natural in the small pitch angle regime, simply 
as a result of decreasing average pitch angle. As 
seen in equation (2) and described in §3.1, when 
decreases, the characteristic SPD frequency v s 
approaches the characteristic IPD frequency v m 
causing the fitted value of a to increase from —2/3 
to 0. Meanwhile, E p oc v m oc sin 1 ? will decrease as 
the pitch angle decreases. We believe this physi- 
cal effect produces the negative correlation seen in 
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Fig. 6. — The E p distribution for those GRBs in 
the IPD regime (solid line) and SPD regime (dot- 
ted line). 

Figure 4 and dominates over the positive instru- 
mental correlation between a and E p . [The in- 
strumental effect is lessened in this regime because 
of the steeper low energy slope (relative to the 
isotropic case) , which allows the spectrum to reach 
its low energy asymptote more quickly.] We point 
out, however, that the quantitative value of the 
slope is somewhat shallower than what is naively 
expected if only the mean value of the pitch an- 
gle ^ changes, while all other parameters remain 
constant. This is illustrated by the dashed line in 
Figure 4, which shows how E p changes as a func- 
tion of a as the pitch angle decreases. To obtain 
this curve, we simulate SPD spectra for a series 
of decreasing pitch angles (given a constant mag- 
netic field, 7 m , and bulk Lorentz factor T), add 
Poisson noise to the spectra, and then fit a Band 
spectrum to each curve; this gives us a value for 
E p and a as a function of decreasing pitch angle. 
Naturally, dispersion in the magnetic field, -f m or 
bulk Lorentz factor will tend to weaken the corre- 
lation (as is observed). We discuss this further in 
§5.1.4 below. 



4.2.3. SAS Regime 

This negative correlation appears to continue 
for the self-absorbed spectra (SAS). The low E p 
value in the SAS regime is consistent with our sug- 
gestion (LPOO) that the self-absorption frequency 
is at the lower edge of the BATSE window (around 
50keV). However, there are few bursts in this 
regime and this may be a result of an observa- 
tional selection effect. For very large values of a, 
bursts with lower values of E p (everything else be- 
ing equal) will have a better chance of triggering 
the BATSE detector. Unfortunately, there is not 
yet sufficient data to test this conjecture. We note 
that such a selection effect does not seem likely 
to explain the negative correlation in the SPD 
regime, because of the small difference in the E p 
distribution between the IPD (a < -0.7) and SPD 
( a <L —0.7) cases, as shown in Figure 6. 

4-2.4- Summary 

The a — E p correlations provide, at least quali- 
tatively, further support for the synchrotron model 
and possibly the small pitch angle distribution 
(SPD) scenario. Our most intriguing result is 
that at a > —0.7, the a — E p correlation goes 
from positive (as expected purely from instru- 
mental effects discussed above and in LPOO) to 
negative. This negative correlation is a natural 
expectation in the SPD emission regime of our 
synchrtron models. However, there are some inter- 
esting quantitative points to address. The dashed 
line in Figure 4 gives an example of how E p should 
change as a function of a in the BATSE spectral 
window, if only the mean of the electron pitch 
angle * changes (all other parameters such as B 
and 7 m remaining constant). The parameter E p - 
in going from the isotropic to anisotropic electron 
distribution regime - should decrease by a factor 
<~ 7m/2, as sin^t decreases from <~ 1 to <~ l/7m 
(see equations 2 and 3). It is often assumed that 
the electron Lorentz factors in GRBs are ~ 100 - 
approximately equal to the bulk Lorentz factors 
which need to be at least this large to keep the 
medium optically thin to pair production and in- 
verse Compton scattering (see Lithwick and Sari, 
2001, for a recent discussion of this issue). In these 
cases, we should see a decrease in E p by a factor 
~ 50 as we transition from the isotropic to very 
small pitch angle regime. Although this strong 
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negative correlation will be weakened to some de- 
gree due to dispersion in the intrinsic values of ^ 
and 7 m , as well as variation in the magnetic field, 
we might expect a stronger decrease in E p than 
what is seen in the data for the large values of j m 
mentioned above. The relatively small observed 
change in E p could mean several things: 

(a) The minimum electron Lorentz factor or the 
magnetic field of the electrons increases as we 
transition to a physical regime in which electrons 
are accelerated primarily along the magnetic field 
lines. This may be a very plausible explanation - 
there may exist physical situations which require 
either a higher magnetic field or characteristic 
electron Lorentz factor, in which it is very effi- 
cient to accelerate along the magnetic field lines. 
The details of this are beyond the scope of this 
paper. 

(b) The electron Lorentz factors could be much 
smaller than 100 - in fact, in an internal shocks 
model, it is expected that the electron Lorentz 
factors will be on the order of the relative Lorentz 
factor of the two shells, which can be on the order 
of a few (Piran, 1999). This would be consistent 
with the relatively small decrease in E p in the 
SPD regime, seen in the data. We caution, how- 
ever, that with electron Lorentz factors so low, we 
require a larger value of the B field (for a given 
observed value of, e.g., E p ) and, in addition, some 
of the assumptions implicit in our spectral models 
(e.g. extreme relativistic velocities) may not hold. 

(c) This model is incorrect and an alternative 
explanation is needed to accommodate bursts 
above the a = —2/3 line of death. For ex- 
ample, the photon indices a between —2/3 and 
the self-absorption value of 1 may simply be due 
to the presence of both the absorption frequency 
and minimum electron frequency present in the 
BATSE window. As long as v m > v a and spec- 
tral fits place the characteristic break energy at 
E p cx v m , then the low energy photon index will be 
a weighted average of —2/3 and 1 (depending on 
the relative values of v m and v a ). We might then 
explain the observed negative correlation between 
a and E p for —2/3 < a < by the following: As 
v m decreases (while v a remains roughly constant 
or even increases), E p decreases and we get less of 
the —2/3 portion of the spectrum relative to the 
slope = 1 portion. This will cause the value of a 
to increase relative to the —2/3 value (of course, 
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Fig. 7. — Photon flux / 7 vs. E p for the binned 
2,026 spectra in our sample. The dotted line shows 
the correlation (for all bursts at the same redshift) 
as a function of an increasing B-ficld. The dashed 
line shows / 7 vs. E p as j m increases. 



in this case, we expect that the Band fits with a 
single break will not be as good). We do not go 
into further detail on this subject and only present 
it as another possibility. 

4.3. Observed Total photon Rux-E p Corre- 
lation 

Figure 7 shows a correlation between the av- 
erage burst total peak flux and E p (we point out 
that this correlation is also present between the 
normalization or height of the spectrum and E p , 
so that the value of the peak flux is not greatly 
affected by the values of the low and high en- 
ergy spectral indices and is more a measure of 
the brightness or overall emission power). The 
correlation is well established and has also been 
reported within individual bursts (Mazets, et al., 
2001). This correlation can be due to several ef- 
fects in a synchrotron model including increasing 
magnetic field, minimum electron Lorentz factor 
j m , and/or bulk Lorentz factor T; in Figure 7, 
we show / 7 as a function of E p in an IPD regime 
for changing magnetic field (dotted line) and min- 
imum electron Lorentz factor (dashed line). We 
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point out again, however, that the theoretical lines 
describe the relation between luminosity and the 
value of E p in the cosmological rest frame of the 
burst. A dispersion in the redshift distribution of 
the GRBs will tend to smear out this theoretical 
correlation. We interpret this correlation as a re- 
lation between E p and flux from pulse to pulse (or 
emission episode to emission episode), and have 
suggested ways to reproduce this by varying the 
plasma parameters in the context of synchrotron 
radiation. We note that such a correlation has 
been reported in the decay phase of individual 
pulses of bursts (e.g., Borgonovo and Ryde, 2000, 
Ryde and Svensson, 2001). Such a correlation in 
a single pulse may arise from a change in plasma 
parameters within a single emission episode, but 
can also be produced by relativistic beaming ef- 
fects alone (photons from the edge of the beaming 
cone arrive later than those from the middle of the 
cone; the edge photons have a smaller doppler fac- 
tor [by 1/r 2 ] and so have both a smaller value of 
flux and E p at later times in the pulse). In any 
case, the above synchrotron interpretation appears 
as a plausible candidate for the description of the 
correlation observed in the data. 

There may be correlations as a result of other 
physical effects, and correlations between the 
physical parameters themselves. For example, if 
particle number were conserved in an emission 
episode, then the normalization of the particle 
distribution would increase as q increased - this 
would cause an increase in the overall normal- 
ization of the spectrum; since the total particle 
number J N(j)dj = J N 1+ ^ / /™]l +q d 1 , then N 
scales roughly linearly with q. 

5. Examples of Spectral Evolution in In- 
dividual Bursts 

Because each pulse is a separate emission 
episode in our model, it is useful to examine what 
our models can tell us about the specific behavior 
of different internal shocks within an individual 
burst. This section is intended primarily to give 
the reader a feel for the different types of spectral 
evolution present in GRBs. All of the following 
spectral fits are again taken from the catalog of 
Preece et al. (1999). In Figures 8-13, we display 
the total photon flux, / 7 (top panels), the peak 
of the vF v spectrum, E p (middle panels), and the 
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Fig. 8. — Evolution of photon flux / 7 (top panel), 
E p (middle panel), and a (lower panel) for burst 
3492. The dashed line in the top panel is the total 
GRB count rate arbitrarily normalized to the peak 
flux value. Note that this burst remains in the IPD 
regime throughout. 

low energy photon index, a (lower panels), as a 
function of time, t, for several GRBs; the dashed 
histogram superimposed on the / 7 (t) plot is the 
time profile (of detector counts) of the burst. We 
now qualitatively discuss the behavior in the con- 
text of our synchrotron emission scenarios: 

1. Burst 3492: Figure 8 displays the spectral 
evolution of burst BATSE trigger number 3492. 
This is a good example of how E p tracks the flux 
(or counts) of the burst. The parameter a some- 
what tracks E p in this case, but there is evidence 
for a change in this trend which could be indica- 
tive of a change in the lower cutoff parameter q of 
the electron spectrum. We note that - according 
to its value of a - this burst remains in the IPD 
regime throughout its duration. 

2. Burst 5567: Figure 9 is another example of 
E p tracking the flux and a tracking E p - entirely 
in the IPD regime. 

3. Burst 2286: In Figure 10, E p appears to 
correlate with the flux, and begins at low or soft 
values and ends at higher, harder values (note that 
the first point in the plot at t = could be indica- 
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Fig. 9. — Same as Figure 8, but for burst 5567. 

tive of a precursor). This suggests that during the 
later emission episodes, either the magnetic field 
or minimum Lorentz factor changed to cause an 
increase in E p . How this is physically achieved 
is a complicated issue, and is related to the level 
of turbulence and/or the specifics of the particle 
acceleration mechansim in the shock. The inter- 
pretation of spectral evolution in our models gives 
us a way of directly interpreting how the physics 
can change from shock episode to shock episode. 

4. Burst 2855: In Figure 11, we see that 
a evolves from about 1, interpreted as a self- 
absorbed situation (and note the low value of E p ), 
to about zero, a value appropriate for the small 
pitch angle case, and then back to 1 or the self- 
absorbed regime at the end of the burst (where 
again E p is at its lowest). If this interpretation is 
correct, the physics required to produce this be- 
havior is intriguing to say the least - particularly 
switching back to the self-absorbed regime at the 
end of the burst. We also point out that for this 
burst - which is entirely in the SPD and SAS emis- 
sion regimes - a appears to roughly anti-correlate 
with E p (see §3.1 and Figure 4 above; although 
admittedly, the error bars are large here and a 
constant E p is a statistically acceptable - although 
not the best - description of the data). 



Fig. 10. — Same as Figure 8, but for burst 2286. In 
the bottom panel, the upper and lower dashed hor- 
izontal lines delineate the SPD and IPD regimes, 
respectively. 

4. Burst 3489: Figure 12 shows an addi- 
tional example in which a appears to roughly 
anti-correlate with E p during the various emission 
episodes of this burst, particularly from t ~ 9s to 
14s. 

5. Burst 1886: Figure 13 shows an example 
of a transitioning from the SAS to SPD regime; 
meanwhile, E p appears to roughly track the flux 
of the burst. 

The different types of behaviors seen in the in- 
dividual bursts are interesting. On one hand, the 
spectral evolution appears to be highly varied and 
all types of behavior can exist (see also Kargatis, 
1994). On the other hand, certain trends appear 
to be present in the data that are at least qualita- 
tively consistent with the global trends presented 
in §4. Of course, we acknowledge that these global 
trends will not be seen in all bursts; establishing 
statistically robust trends for each particular emis- 
sion regime in each particular burst in an attempt 
to say something about the average trends among 
all bursts is not statistically meaningful (see sec- 
tion 4 for discussion of this issue). We do em- 
phasize again that all types of behavior can occur 
within individual bursts, and wc can use our mod- 
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Fig. 11. — Same as Figure 8, for burst 2855. The 
horizontal dashed line separates the SAS (above 
the line) and SPD (below the line) regimes. Note 
that this burst remains entirely in the SPD and 
SAS regimes througout its entire duration. There 
is also a rough anti-correlation between a and E p 
for this burst. 

els to interpret the physics behind any particular 
evolution. For example, Crider et al. (1997) find 
that in GRB910927 6 , there is a positive correla- 
tion between a and E p in the SAS regime (and 
very slightly in the SPD regime, although there 
are only 6 points here) . This burst appears to have 
one (or at most, two) smooth, broad pulses, so 
that there is probably only one (or two) emission 
episode(s) giving rise to the various spectra. A 
positive correlation in this case, could be due to a 
number of factors in our model. Of course, as dis- 
cussed in LP00 and in §2, the finite bandwidth of 
the BATSE spectral window will tend to produce a 
positive correlation between a and E p . However, a 
physcial explanation is also plausible; for example, 
it is possible that the the self-absorbed (or small 
pitch angle) portion of the spectrum is decreasing 
relative to the IPD (a = —2/3) portion due to a 
decrease in, say, the magnetic field in the shock. If 

6 Note that this burst is not present in the Prcece et al. cat- 
alog because it does not meet one of the catalog's selection 
criteria - namely, a peak flux greater than 10 ph/cm /s on 
the 1024ms timescalc. 



Fig. 12. — Same as Figure 8, but for burst 3489. 
During the various emission episodes (to the ex- 
tent they can be delineated), E p is on average 
lower for the higher values of a in the SPD regime. 

E p corresponded to the self-absorption frequency, 
this would produce a positive correlation between 
E p and a. This correlation would be exacerbated 
if the cutoff to the low end of the electron distri- 
bution evolved to a more flat distribution (lower q 
values) throughout the pulse. 

Our point is that no matter what the behavior, 
our models can be used as diagnostics to interpret 
what types of physical changes may be occuring 
from pulse to pulse. Note that if we can in fact 
characterize our synchrotron emission regime by 
the value of a, then it is particularly interesting 
when we see bursts switch emission regimes from 
pulse to pulse, (as a caveat, however, we note that 
the error bars are significant on some of the spec- 
tral fits and a clear cut interpretation is not always 
readily available). This suggests that the funda- 
mental plasma physics (i.e. the particle acceler- 
ation) can vary depending on the internal shock 
conditions within a single burst. 

6. Discussion and Conclusions 

In this paper, we have explored the validity of 
the synchrotron model in explaining the behavior 
of the time resolved GRB spectra. In our model, 
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Fig. 13. — Same as Figure 8, but for burst 1886. 
The horizontal dashed line delineates the SPD 
emission regime. This burst evolves from an SAS 
to SPD regime. 

there are three different emission regimes, all dis- 
tinguished by the value of the low energy photon 
index a. Our model also accounts for the instru- 
mental correlation between a and E p produced 
when E p is near the lower end of the B ATSE spec- 
tral window (in this case, the spectrum has not 
reached its low energy asymptotic value and so a 
will be softer or lower for smaller E p 's). In the 
IPD regime, the synchrotron emission is optically 
thin from an isotropic electron distribution, and 
the low energy photon index a is less than or equal 
to —2/3. In the SPD regime, the emission is pri- 
marily from electrons with small pitch angles and 
the spectrum consequently has a value of a ~ 0. 
Finally, in the SAS regime, the synchrotron pho- 
tons are self-absorbed and there is a steep cutoff 
(a = 1 or 3/2 depending on the relative values of 
the absorption and minimum electron frequency) 
in the low energy spectrum. 

We have shown that these models provide an 
excellent description of the existing data, and 
have presented spectral fits for a few bursts. Our 
results suggest that the Band parameter a is a 
good diagnostic of the relevant emission regime. 
We also have presented the spectral evolution of 



the Band spectral parameters (specifically photon 
flux, E p , and a), and interpret the behaviors in 
the context of our model. We point out that the 
behavior is quite varied and there is often evi- 
dence that a burst switches emission regimes from 
pulse to pulse. In attempt to characterize general 
evolutionary trends exhibited in GRBs, we have 
combined all 2,026 separate time resolved spectra 
(from the 80 bursts in our sample) and looked for 
correlations in this data. Our main results are as 
follows: 

• We find that the majority of the bursts lie in 
the IPD regime (~ 55%), but that a significant 
fraction of the bursts (~ 40%) are in the SPD 
regime. This has interesting implications for the 
particle acceleration studies in GRBs as we discuss 
below. We find that only a small fraction (~ 5%) 
of bursts are in the self-absorbed (SAS) regime. 

• We show that there is a strong correlation be- 
tween a and E p in the IPD regime, which can be 
attributed to the positive instrumental correlation 
discussed in LP00. However, there appears to be 
a negative correlation between a and E p in the 
SPD regime (albeit with 4cr significance). We in- 
terpret this as evidence of the effects of decreasing 
electron pitch angle as a transitions from —2/3 
to 0. This physical effect overwhelms the positive 
instrumental correlation that is dominant in the 
IPD regime. However, for this physical effect to 
accommodate the data quantitatively, we need ei- 
ther an increase in the minimum electron Lorcntz 
factor as we transition to the SPD regime, or for 
the absolute value of the electron Lorentz factor to 
be only on the order of a few. This may very well 
be the case in an internal shocks model, where 
the relative Lorentz factor of the two shocks is the 
relevant scale for the electron Lorentz factors. 

• We find a strong positive correlation between 
the photon flux and E p , which can be explained 
by several effects in a synchrotron emission model. 
Most notably, a change in the magnetic field or 
minimum Lorentz factor will produce this type of 
a correlation. 

Our results bring to light the fact that parti- 
cle acceleration in GRBs is a quite poorly under- 
stood problem. Usually, it is assumed that the 
radiating particles in GRBs are accelerated via 
repeated scatterings across the (internal) shocks. 
This is because shocks can quickly accelerate par- 
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tides to very high energies, through repeated scat- 
terings across the shock (the scattering agent be- 
ing plasma turbulence from, e.g., a two-stream in- 
stability Medvedev and Loeb, 1999). This mech- 
anism, however, predicts several features in the 
electron distribution not borne out by the data. 
First, it has been shown (Kirk et al., 2000 and ref- 
erences therein) that these repeated crossings of 
the shock result in a power law particle distribu- 
tion with a well defined index, p = —2.23, which 
would give a high energy synchrotron photon in- 
dex (3 of -1.62 (or -2.12 for the "cooling" spectrum, 
e.g. Sari et al., 1997). Although this is consistent 
with some afterglows, this is certainly is not true 
for many bursts in the prompt phase. In our syn- 
chrotron models above, the high energy photon 
index (3 = — (p + l)/2, where p is the high energy 
index of the emitting particle distribution. The 
parameter (3 can vary by a factor of 4 (or more!) 
throughout a single burst (see, e.g., Preece et al., 
1999), reflecting a huge variation (from 1 to 9) in 
the parameter p of the underlying particle distri- 
bution - this is well beyond the statistical limits 
placed on p by shock acceleration simulations. In 
addition, shock acceleration predicts an isotropic 
distribution of electrons. Our work suggests that 
in a large fraction of GRBs, the particle accelera- 
tion is not isotropic but along the magnetic field 
lines. 

In fact, when the Alfven (phase) velocity (3a = 
B/^nnmpC 2 ) 1 / 2 (in units of the speed of light) 
is greater than 1 , stochastic acceleration becomes 
more efficient than shock acceleration (Dung and 
Petrosian, 1994). It is well known that in this 
case the electric field fluctuations 5E ~ (3a5B ex- 
ceed the magnetic field fluctuations, which means 
a faster acceleration (oc SE 2 ) than scattering rate 
(oc SB 2 ). Under these conditions, once the parti- 
cle crosses the shock front into the turbulent region 
behind the shock, it will undergo stochastic accel- 
eration much faster than it can be turned around 
to cross the shock again. This is the situation 
for GRBs which have inferred magnetic fields of 
B <~ 10 5 G, and densities of n < 10 8 cm~ 3 so that 
(3a 1- The shape of the spectrum is then de- 
termined by the relative values of the diffusion co- 
efficients and the rates of energy and pitch angle 
changes resulting from the interaction of the in- 
jected particles with the plasma turbulence behind 
the shock. Hence, we emphasize the importance 



of an investigation of particle acceleration in GRB 
internal shocks. 

Lastly, we would like to emphasize the impor- 
tant role upcoming (and current) GRB missions 
will have in laying to rest the questions raised in 
our investigation. With the launch of HETE-2, 
and the upcoming launches of Swift and GLAST, 
we will be able to get high quality broadband 
(from a few eV to GeV) spectra of the prompt 
GRB emission; such spectra will allow us to test 
our models more stringently and constrain all of 
the emission mechanisms that may play a role in 
GRBs (for example, we may see a synchrotron sclf- 
Compton component in the GLAST energy range; 
Dermer et al., 2000). Data from these satellites 
are sure to shed significant light on the photon 
spectrum and particle acceleration mechanisms in 
Gamma-Ray Bursts. 
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